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' We present the results of the tree-level calculation of deep-inelastic leptoproduction, in- 

, eluding polarization of target hadron and produced hadron. We also discuss the dependence 

■ on transverse momenta of the quarks, which leads to azimuthal asymmetries for the produced 

' hadrons. 
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cn ; 

' In recent years several possible ways to probe the structure of hadrons in hard scattering processes have been 
ly-j , pointed out. Most well-known are the inclusive lepton- hadron [IH) scattering experiments that provide detailed 
^\ information on the unpolarized and polarized quark distributions. This information is valuable as it, within the 
framework of Quantum Chromodynamics (QCD), can be expressed as well-defined matrix elements of quark 
and gluon operators within the nuclcon. For instance, the operator product expansion relates moments of quark 
1 distributions to expectation values of local operators that represent static properties of the nucleon such as 
Oh, its electric or axial charge, baryon number, etc. Some of these matrix elements can also be determined in an 
^ independent way, such as the axial charge of the nucleon from neutron decay. We will consider in this paper 
matrix elements of nonlocal combinations of quark fields . They represent forward quark-target scattering 
amplitudes that can be interpreted as quark momentum distributions or multi-parton distributions . 

In inclusive ^H scattering one can determine combinations of the quark distributions fi{x) and gi{x) where x 
is the fractional lightcone momentum x = p'^ / of a quark with momentum p in a hadron with momentum P. 
d The index a refers to the flavor of the quark. The antiquark distributions found in inclusive £H scattering are 
denoted fi{x) or f1{x), etc. These quark distributions appear in the structure functions, which in the Bjorken 
limit scale in the variable x^ = /2P ■ q, where = —q^ with q the (large) momentum of the virtual photon. 
Restricting ourselves to the exchange of photons, the experimentally accessible leading structure functions are 

E^k^^ ^ 2F,{x.,Q-) if^ixs) + m^u)) , (1) 

a 

2g^{x,,Q^)^Y.^l (2) 

a 

where Ca are the quark charges. The functions appear in leading order in l/Q and are referred to as twist- 
two functions since their moments are matrix elements of local operators of twist two. In naming the quark 
distributions we extend the scheme proposed by Jaffe and Ji |^,^. For the unpolarized quark distributions, 
often denoted as /"'{x) or a{x), the notation becomes fi{x). For the polarized quark distributions, for which 
one finds A/''(a;) or Aa{x), the notation becomes gi{x), 

The only other distribution function that can be obtained in inclusive deep-inelastic IH scattering up to order 
l/Q is the twist-three function g^{x). It appears in the structure function 

2 g^{x^, Q2) (ff2 (^b) + 92{^b)) , (3) 

a 

where (^0 = ~9i{x) + gxi^)- The functions can be expressed as matrix elements of quark and gluon field 
operators |o,0| and their moments involve local operators of twist three M . 
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In processes involving at least two hadrons, it is possible to measure one more twist- two quark distribution 
function. This function, the transverse spin distribution hi(x) (also known as ATa{x)) and the antiquark 
distribution hi (x) are chirally odd, hence they must be combined with some other chirally odd structure PpO) . 
This is possible in Drell-Yan scattering or in semi-inclusive £H scattering. In the latter case the chirally odd 
structure is found in the quark hadron fragmentation process. At the twist-three level, two more chirally 
odd quark distributions, e°'{x) and h1{x) appear 

In processes involving at least two hadrons or two jets, one also has the possibility to investigate azimuthal 
dependences. This involves a momentum that is transverse to the large momentum flow, e.g., the azimuthal 
dependence of hadrons around the (spacelike) direction defined by the virtual photon. This is the specific 
case investigated in this paper. Other examples are the azimuthal dependence of the produced muon pair in 
Drell-Yan scattering or the azimuthal dependence of hadrons around the jet-jet axis in two-jet events in e^e~ 
annihilation. 

As has been discussed first by Ralston and Soper |Tl[] and recently by us , one has the possibility to measure 
in total six quark distribution functions depending on the lightcone fraction x and the transverse momentum 
squared p^, which depends on x and the quark virtuality p^. In this paper we investigate the full tree-level result 
for semi- inclusive £H scattering up to order 1/Q which in addition to the six twist-two distribution functions 
involves eight twist-three functions. After pr-hrtegration only the three twist-two and three twist-three {x- 
dependent) functions discussed above remain. For semi-inclusive scattering one also needs to study the quark 
fragmentation process involving a number of twist- two and twist-three fragmentation functions. The treatment 
of the fragmentation functions proceeds along the same lines as the distribution functions with the exception 
that one also can have some time-reversal odd structure. This structure emerges as the outgoing hadron in £H 
scattering is not a plane wave state in contrast to the incoming hadron. 

An important motivation to study quark and gluon distribution and fragmentation functions in hard scattering 
processes are the existence of factorization theorems in QCD that imply the universality of these functions, 
notwithstanding their (logarithmic) scale dependence. Such factorization theorems exist for the p^-iiitegrated 
functions |^^. They do not exist for the unintegrated distributions [ p^ljT^ ]. Nevertheless, we consider the present 
complete tree-level result up to (subleading) order 1/Q as an important step. It represents the extension of the 
naive parton model result to subleading order and shows the dominant structures to be expected in the cross 
section although QCD corrections, such as Sudakov effects, may affect the asymmetries iQ. At the tree-level 
the only QCD input at order l/Q is the use of the equations of motion which assures the electromagnetic gauge 
invariance. In view of possible future measurements with electrons or muons at DESY, CERN or Fermilab 
it is important to be aware of possible azimuthal asymmetries, not only because they contain valuable new 
information on the quark and gluon structure of hadrons but also in order to estimate corrections to Pt- 
integrated quark distributions coming from variations in angular acceptance. 

The structure of this paper is as follows. In section 2 we present the formalism of £H scattering followed in 
section 3 by the analysis of quark distribution and fragmentation functions. This is followed by the full tree- level 
result up to order 1/Q in section 4, the basic result of this paper. The observable consequences of this result 
are outlined in the next two sections, starting with the full py-integrated result for £H — > hX in section 5. The 
cross sections are given including polarization of lepton and target hadron and including measured final state 
polarization, all of this for the case of spin 1/2 particles. In section 6, the result for azimuthal asymmetries 
involving unintegrated distribution and fragmentation functions is presented, including only polarizations of 
initial state particles. Finally the results are summarized in section 7. 

II. FORMALISM 

We consider deep-inelastic semi- inclusive lepton-hadron scattering, £ + H £' + h + X, in which a lepton with 
momentum / scatters off a hadron H with momentum P and one hadron h with momentum Ph belonging to 
the current jet is measured in coincidence with the scattered lepton with momentum I' . The momentum of the 
exchanged virtual photon q — I — I' is spacelike with —q^ = Q^ ^ oo. We will use the invariants Xb = Q^/2P-q 
~ —Ph ■ q/Ph ■ P, y ^ P ■ q/ P ■ I ^ Q^/xbS (where s is the invariant mass squared of the photon-hadron 
system), and Zh = P ■ P^/P ■ q w —2Ph ■ q/Q^, the approximate sign indicating relations valid up to order 
1/Q^ corrections. Furthermore, we use a transverse (two-component) vector which is small (i.e., of order 
M). In the frame where h and H are collinear it is precisely the transverse component of q, while in the frame 
where P and q are collinear it is up to 1/Q^ corrections equal to « —Ph^/zh. Note that we will henceforth 
systematically neglect 1/Q^ corrections even if we use equal signs. 

The cross section for semi-inclusive lepton-hadron scattering is given by 



^^;^^^^ = i|...i^..2M>V- (4) 

The leptonic tensor is standard and is for polarized leptons (neglecting the lepton masses) given by 

L^^IX; I'X') = 6xx' i^yl + 21J'^ - Q^g^, + 2iX e^^p^q^r) , (5) 
with the lepton helicity being A = ±1. The hadronic tensor is given by 

2MW^.(?; PS; PhSh) = j^T. J (2^)3^po + P - Px - Ph) 

X {PS\J^mPx;PhSh){Px;PhSh\MO)\PS). (6) 

The role of the polarization vectors 5* and Sh is discussed in Appendix A. For the inclusive process one finds 

da -K o? 
dxe dy Q"* 

with the hadronic tensor given by 

X 



yL^,2MW'\ (7) 



2MW^.{q;PS) = ^Y. J (2tf2P^^ + P- Px){PS\MO)\Px){Px\MO)\PS) 

= ^ y d'x e^'-'' (PSliMx), MO)]\PS}. (8) 

In order to expand the leptonic and hadronic tensors it is convenient to work with vectors orthogonal to q. 
A spacelike vector is defined by and a timelike vector defined by P'^ = P'^ — {P ■ q/q^) q^, 

= (9) 

= — P** = 2xb P^ + q", (10) 

P • q 

These vectors satisfy T'^ = Q'^, Z'^ = -Q'^. We will consider normalized vectors = T^'/Q, z^" = Z^'/Q = 
—q'^ = —q^/Q. In the space orthogonal to z and t one has the tensors 

g^^- = gf^- + ^f^^- -i'^i'', (11) 

ef = ei^'^'H^q, = ^ e'^^^'-P^?,. (12) 



One vector in the perpendicular space is Pj^j^ = g'"' Phv We will frequently use the normalized vector /i'' = 

pi:j\phi^\- 

Also the lepton momenta / and V = I — q can be expanded in t, z and a perpendicular component, 



y 2 2 2y y 

where f = is a spacelike unit-vector in the perpendicular direction lying in the (lepton) scattering 

plane. The kinematics in the frame where virtual photon and target are coUinear (including target rest frame) 
is illustrated in Fig. 1. With the definition of x, we have for the leptonic tensor 



r IJtV ^ 



-2(l-y + iy')5r + 4(l-2/)t>t^ 



+4(1 - y) (^x^x'^ + i ^ + 2(2 - y) t^^'x-'^ 



-iXy{2 - y) e^^ - 2iXy^fY^y {["e^^^'f 



(14) 



Theoretically, it is often convenient to work in the frame in which the two hadrons H and h are eoUinear. 
Except for corrections proportional to M"^ /Q"^ the vectors /Q and P{^/Q determine two lighthke directions 
and n'^, respectively. Using for momenta the representation p — [p~ ,p'^ ^p^], where p^ = [p^ ±p^)/'\/2, one 
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obtains in a frame in which H and h are collinear 
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(15) 
(16) 
(17) 



where A fixes the frame and n 

momentum q-r with length = — = — Q^. It is easy to see that P^'^j^ 
Transverse projection operators are given by 



n_ = 1. In a collinear frame as given here, the photon acquires a transverse 
2 _ _n2 Tf ^c, .00,. PA' - and thus = \PhA/zh. 



n_i n , 



(18) 
(19) 



where the brackets around the indices indicate symmetrization of these indices. Note that these transverse 
projectors are not identical to the perpendicular ones defined above if the transverse momentum of the outgoing 
hadron does not vanish. The lightlike directions, however, can easily be expressed in t, z and a perpendicular 
component. 



^ [t'' + z^- 

1 
71 



Q 



(20) 
(21) 



showing that the differences are of order 1/Q. Especially for the treatment of azimuthal asymmetries, it is 
important to keep track of these differences. 



III. QUARK CORRELATION FUNCTIONS 



In this section we investigate the 'soft' hadronic matrix elements that appear in the diagrammatic expansion 
of a hard scattering amplitude |l7j . The most important matrix elements that can be measured in deep-inelastic 
semi-inclusive lepton-hadron scattering are of two kinds. The first is 



'^viP^S-^p) 



(27r)4 



d^x e^P-'^iP, S\'ipj{0)CiO,x; path)iPi{x)\P, S), 



(22) 



where a summation over color indices is implicit and flavor indices are omitted. This expression, referred to as 
a quark-quark correlation function, is a forward matrix elements of quark fields with explicit Dirac indices i and 
j in a hadronic state characterized by its momentum P and spin vector S. The quark fields are accompanied 
by a path ordered exponential (link operator). 



£(0, x; path) = V exp 



ds^^A^is] 



(23) 



necessary to render the definition color gauge- invariant. It introduces, however, a path-dependence. In order to 
obtain a simple parton interpretation (Fig. 2) for the correlation functions the choice of path is important, as 
we will see later. Technical details will be discussed in Appendix B. The other type of matrix elements measured 
in semi-inclusive lepton-hadron scattering are the quark decay functions [D (Fig. 2) 



^^3{Ph,Sh■,k)^Y.Ji^ J e""" {0\mx;P^th)Mx)\Ph,Sh;X){Ph,Sh;X\ijj{0)\0) 



d'^x e**--^ (0|/:(0,a;;path)V'z(a;)a;a,,Vi(0)|0), 



(24) 



where an averaging over color indices is implicit. Again a quark link operator accompanies the quark fields. 
The second way of writing using hadron creation and annihilation operators is to be considered as a shorthand 
of the first expression, involving (out-)states \Ph, Sh', X) where Ph and Sh are the momentum and spin vectors 
for the produced hadron. 

The correlation functions are constrained by the hermiticity properties of the fields, and by parity and time 
reversal invariance. Without constraints from time reversal invariance one has 



(25) 



+ ^7 [v ■ S) + A^{p- S) + i Ag a^^^jz S^P, + i Aw a^^js S^p, 
+1 An {p ■ S) p^Pu + ^12 V^", 



where the first four terms do not involve the hadron polarization vector. Hermiticity requires all the amplitudes 
Ai = Ai(p ■ P,p^) to be real. The amplitudes A4, A^ and A12 vanish when also time reversal invariance applies. 
A similar decomposition in 12 amplitudes can be made for A, but because of the fact that out-states appear in 
the definition of the correlation functions time reversal invariance cannot be used and none of the amplitudes 
vanish. Because the time-reversal condition for A4, A^ and A12 reads A* ~ ~Ai we refer to these amplitudes 
as time-reversal odd. 

The above matrix elements are assumed to be nonzero only for quark momenta limited to a characteristic 
hadronic scale, which is much smaller than Q^. This means that in the above matrix elements p^, k'^, p ■ P 
and k ■ Ph <C Q^. With the definition of momenta in the coUinear frame (Eqs [l^ - one then arrives at the 
Sudakov decomposition, for the quark momenta. 



p : 

k - 



2x P+ 



2p: 



xQ 
Xb\/2 

ZhQ 

zV2 



Z(fc2+fc2) 



Pt 



Zh 



Q\/2 



1 

-Ph 

z 



(26) 
(27) 



where in the last expressions for p (k) the terms proportional to n_ (n+), irrelevant in leading, 0(1), or 
subleading, 0{1/Q), order have been omitted. It turns out that up to subleading order only the combinations 
J dp~ ^(P, S\p) and J dfc"*" A(P;i, Sh\ k) are important. For cross sections integrated over transverse momenta 
only the combinations J dp^cfpT $(P, S;p) and J dk^cPk^ ^{Ph, Sh; k) are of importance at leading order. 

The functions that appear in hard scattering processes can be expressed as specific Dirac projections of the 
correlation fimctions, integrated over p~ , 



$Pl(a;,p^) =\j dp- Tr($r) 



2 (27r)3 



e'P-« (P,5|^(0)r/:(0,C;n_)V(OI^,5> 



(28) 



C+=o 



In general the matrix element contains an infinite tower of operators because of the link operator. Choosing the 
path in the definition of <I> to lie in the plane = and in essence along the direction x~, indicated with 'path' 
— ri- (see Appendix B), one finds that in an appropriately chosen gauge in leading order only one (bilocal) 
operator combination contributes, which allows a parton interpretation. 

The projections ^['"l depend on the fractional momentum x = p~^ / P^ and on p^ and furthermore on the 
hadron momentum P (in essence only P+ and M). Depending on the Lorentz structure of the Dirac matrix T 
the projections $['"1 can be ordered according to powers of MjP^ multiplied with a function depending only on 
X and p\. = —p\- It is easy to convince oneself that each factor MjP^ leads to a suppression with a power MjQ 
in cross sections iQ. In analogy with inclusive scattering we therefore refer to the projections as having a 'twist' 
t related to the power {AI / P^)*--"^ that appears. With this definition the moments (in x) of the p^-integrated 
functions indeed involve local operators of twist t 

It is convenient to use for the spin vector S also a Sudakov decomposition. 



\M XP^ 



2P+ ' M 



XQ 



MXnV2 



XMxj= 
QV2 



(29) 



with X^ + S^ ~ 1 and again neglecting the contribution (proportional to n_) which is irrelevant up to subleading 
order. The following distribution functions appear at leading order in (A//P+)° (twist-two), where as said before 
we omit fiavor indices. 



$[^+l(a;,p,)=/i(x,p^), 

= XgiL{x,pl)+giTix,pl) 



{Pt ■ St) 

M 



XhiL{x,pl) + hiT{x,pl) 



{Pt ■ St) 



M 



Ad* 

= S't hi{x,pl) + -j^ hiL{x,pl) + 



htT{x,pl), 



M2 



with hi = hiT + (p?./2M^) /i^'-j,. At subleading order {M/P+)'^ (twist-three) one finds 
$W(x,pr) = ^e{x,pl) 
^^''\x,pT) = ^f^{x,pl), 
*[^'-l(^,P.) = ^5^(^,P^) + ^ 



^9l{x,pI) +9t{x,pI) 



{Pt ■ St) 

M 



MS: 



^ gT{x,pl) + ^gi{x,pl) + 



S^p't — PtS^ _|_ 



(pipL - hpl^ij) s't 

9t{x,Pt), 



M2 



P+ 



h^{x,pi) 



M 
P+ 



XhL{x,pl) + hT{x,pl)^'^'^ ^'^^ 



M 



with grp = g'j, + (p^/2M^) g^. Prom here on we will often use the shorthand notation 

/ \ _ \ / 2 \ , (Pt ■ St) , 2 \ 

g\s\x,pT) = XgiL{x,PT) H giT{x,PT), 



(30) 
(31) 



(32) 

(33) 
(34) 

(35) 
(36) 
(37) 

(38) 



and similarly shorthand notations h^^, gj- and hg. The results after integration over Pt are lightcone correlation 
ftmctions $Pl(x) depending only on one lightcone momentum component p+ = xP+, for which the nonlocality 
is restricted to x^ , i.e., a;"*" = and xt = Ot- Only the functions /i, gi = giL, hi, e, gx and hL are nonvanishing 
upon integration over Pt • Upon further integration over x ti locs-l matrix element is found, e.g., relating the 
valence distributions to the quark numbers, / dx (/f (x) - /f (x)) = (P, S'|V'"(0)7+'(/'"(0)|P, S')/2P+ = ria - 

In naming the functions, we have extended the scheme introduced by Jaffe and Ji. All functions obtained 
after tracing with a scalar (1) or pscudoscalar (jor,) Dirac matrix arc given the name e.., those traced with a 
vector matrix (7'^) are given the name /.., those traced with an axial vector matrix (7'^75) are given the name 
g,, and finally those traced with the second rank tensor (ia^'^'y^) are given the name h,,. A subscript '1' is given 
to the twist-two functions, subscripts 'L' or T' refer to the connection with the hadron spin being longitudinal 
or transverse and a superscript ' _L' signals the explicit presence of transverse momenta with a noncontracted 
index. 

The twist-two distribution functions have a natural interpretation as parton densities. Using the good and 
bad components of the quark fields, 'iIj± = P±xj} with P± =7^7+72 one finds that the operator combination in 
$[7+] 

is just a density, The other twist-two correlations functions involve in addition chiral projections 

of the quark field, ^r/l = Pr/liJ' with Pr/l = (1 ± 75)/2 or transverse spin projections ip^^^ = P^/^tp with 



T/i 



(1 ± 7^75)/2. The correlation functions then become. 



Jl ' 



.(T) (i) 
/l + /l > 



(J) [7+75] _ j(-R) 



(39) 
(40) 
(41) 



where f[^^ is a matrix element involving the fields V'+rV'+h) etc. Thus, the correlation function $I'''^](x,Pr) 
is just the unpolarized quark distribution, which integrated over Pt gives the familiar lightcone momentum 
distribution $1''' l(x) = /i(x). The correlation function '^^\x,Pt) is the chirality distribution (for massless 



quarks helicity distribution). It can only be measured in a polarized target, as is evident from the presence 
of A and St in Eq. ^ In a transversely polarized (spin 1/2) target the chirality distribution is proportional 
to the transverse momentum of the quark along St- Integrated over Pt only one distribution survives in the 
correlation function, $[''^'''^1 (x) = \gi{x), which can only be measured in a longitudinally polarized target. 
The correlation function <I>[*'^ '''^l(x,pr) is the transverse spin distribution and can also only be measured in 
a polarized target. The possible dependence on target spin and transverse momenta involves three functions, 
as can be seen in Eq. |3^. Using the second expression one sees that one term (involving hi) is insensitive to 
the direction of the quark transverse momentum; this part of the correlation function can only be measured in 
a transversely polarized target. The other two terms in Eq. ^ depend on rank 1 and rank 2 combinations of 
the transverse momentum of the quark and are accessible in longitudinally and transversely polarized targets, 
respectively. Integrated over Pt only one transverse spin distribution survives, 

$[»^'+75](a,) = S'j,hi{x), which 

can only be measured in a transversely polarized target. It involves a chirally odd Q| operator structure, which 
in a cross section can only appear in combination with another chirally odd operator such as a quark mass 
term in inclusive scattering, a chirally odd fragmentation function in semi-inclusive scattering or a chirally odd 
antiquark distribution in Drell-Yan scattering. 

The twist-three distributions cannot be expressed as densities or differences of densities. They are related to 
matrix elements containing one good and one bad quark field. The bad field can in principle be replaced by 
a good field and a transverse gluon field using the equations of motion. Up to 0{1/Q) one anyway needs to 
include quark-quark-gluon matrix elements. For this one needs to consider 

$S,,.(P,5;p) = ^ j d^x e^^-%P,S[4,^{Q)C{Q,x-n_)iD'^{x)U^)\P,S). (42) 

$^^^.(P,5;p) = ^ j d^x e'f-(P,5|^,(0)zD"(0)/:(0,x;n_)i^,(x)|P,5), (43) 

where — id" + gA'^ . The twist analysis is performed by considering the projections 

^f\x,pT)=\ j dp- Tr($Sr) 



^^0^ e^^-« (P, 5|^(0)r £(0, ^; n_ ) iD" [OmW. S) 



(44) 



{+ = 



where C is the same link operator as the one used for the quark-quark correlation functions. Because of the 
choice of link, which always lies along the minus direction (except for the points ^~ = ±oo), one has for the 
correlation function with a longitudinal gluon field the relation 



2{2n) 



p, 5|^,(o)r/:(o, e; n_) iD+MO\P, s) 



= p+$[^l(.T,pr). (45) 

C+=o 



The projections obtained for the quark-quark-gluon correlation functions with transverse gluon fields are not 
all independent from the ones defined for the quark-quark correlation functions, either. Some can be connected 
to quark-quark correlation functions with one good and one bad quark field using the QCD equation of motion, 
[ip — m)ijj{x) = 0. This gives the relations 

gTaf^'i''^}"'*^ ^eTaf^^f"'^^'^ =^iMxe~mh)+eT^JP^S'rxh^^ (46) 
grafi $^^'"'^^=1 = Mx hs - mgu, (47) 
gf - lef = :,f^ _ »e? V/3 (^5.^ - ^ h^) - lefSrp {Mx g'^ - m H^t) ■ (48) 

A useful quantity to consider is the (color gauge-invariant) correlation function 

<^f\x,pT) ^ ^^[^'(x,^^ -p"$F1(x,pO. (49) 

From Eq. ^ one sees that 'I'J''"' = 0, while for the transverse indices a = 1 or 2 it reduces after gauge fixing 
to the pure quark-quark-gluon matrix element (see Fig. 3), 

' , (50) 



^f\x,PT)^g j '^^j^ e^''-Hp,s\m^ A'^{om\p,s) 



needed to calculate the leptoproduction cross sections in subleading order. We use the correlation function 
<I>^ to identify interaction- dependent combinations in the distribution functions. As an example, from the first 
relation above (Eq. we obtain 

^toP^'a^'' ''^^ — i {Mx e — m fi) ~ Ej-ij p^S^ (Hit — xh^) = i Mx e + Sj-ij p^S^ xh^. (51) 

In this way we can rewrite all twist-three functions in a part containing twist-two distribution functions and 
an interaction-dependent part, e.g., e = (m/Mx) fi -\- e. The results for all twist-three functions are explicitly 
given in Appendix C, including results for the p^-integrated functions. 

The antiquark distribution functions in a hadron are obtained from the matrix elements 

'^^,{P,S■,p) = J d^x e-'P--'{P,S\C{0,x;n.)M^)i;j{0)\P,S). (52) 

The antiquark distributions should be defined consistent with the replacement ip ^ Tp'^ — Cip^ , or ^^^^ = +$^['"1 

for r — 7^, i(7fj,uj5 and ^75 and $ = —$'^['"1 for F = 1 and 7^175. Finally, the anticommutation relations for 
fermions can be used to obtain, within the (connected) matrix elements, the symmetry relation 

$„ (P, S; k) = -$y(P, S; ~k). (53) 

For the distribution functions this gives the symmetry relations /i(x,p^) = — /i(— a;,Pr) and identically for 
3iT, hiT, h^Ti 9l and ft-L, while giL{x,pl) = 9il{-x,pI) and identically for h^^, e, /-■-, g'j,, g^, and hx- 

The procedure for the fragmentation functions obtained from the matrix element A(P/j, Sh\ k) is similar. One 
again starts with the general decomposition in twelve amplitudes Ai{k'^,k ■ Ph), replacing P Pt and S ^ Sh- 
in this case time reversal invariance cannot be applied, as the matrix elements A^j involve out-states containing 
a hadron h, in contrast to the plane wave states in the case of the matrix elements $y . A consequence is that 
there are more fragmentation functions. The functions that appear in hard scattering processes can again be 
expressed as specific Dirac projections of the correlation functions, integrated over fc"*". 



A[^l(z,fc;) ^ dk+ Tr(AF) 

d^+d^£,T 



k-=p-/z, kT 



4z(27r)3 



e'^< Tr(0|/:(0, C; n+) ^j{i)alahmm 



(54) 



The arguments of the functions are the lightcone fraction z — Pf^ /k^ and the transverse momentum A;^ — —zk^, 
which is the perpendicular momentum of the hadron h with respect to the quark momentum. In the definition of 
the correlation functions A the path is chosen in the plane x~ = 0, in essence along the x^-direction, indicated 
with 'path' = n-(- (see Appendix B). The spin vector is now parametrized as 



Sh = 



— — j=n^ + n+ + ShT ~ A/t — — + ShT- (55) 



In the last expression the contribution proportional to 71+ , which is irrelevant up to subleading order is omitted. 
The role of the final state spin vector in determining the polarization of the produced hadrons is explained in 
Appendix A. The following fragmentation functions appear at leading order (M/j/P,J^)° (twist-two), 

A['^-l(z, K) = D,{z, fc?) + '-^M^f^ Dixiz, k% (56) 

A[^"^=l(z,fc;) = Gi,(z,fc;) (57) 

A[-'-^=l(z,fc;) = SyHMz.K') + ^Ht{z.k',) +'^Ht{z,k'^). (58) 

Mh Mh 

At subleading order {Mh/P[^Y (twist-three) one finds 

AW(z,fc^) = ^ii;(z,fc;2), (59) 

^h 



A['>''l(z,fc;) - -^Z?^(z,fc;2), (60) 



A\^^^\z,K)^^Es{z,K), (61) 
(z, fc;) = G^(z, k'^) + ^ G^(z, fc;), (62) 

A[-"+^=l(z,fc;) = ^iJ,(z,fc;). (64) 
In the above expressions we have used the shorthand notations Gis etc. standing for 

Gi,(z, fc^) ^ A;, G,l{z. fcg) + , f '"^^ Git(^, k'^). (65) 

There are more functions allowed in this case because of the non-applicability of time-reversal invariance, allow- 
ing more amplitudes in the expansion for Ay as compared to that for <I>ij (see discussion following Eq. ^5|) . The 
fragmentation functions which have no equivalent in the distribution functions, -Dj*^, , El, Et and H , will be 
referred to as time-reversal odd. The results after integration over k'-^ define the lightcone correlation functions 
AP1(z). Nonvanishing upon fc'^-integration are the functions Gi ~ Gil, Hi = Hit + (^t/^-^^/?) ^it' 
El,Gt = G'j, + {kl/2Ml) G^, H and Hl and D^^^^ = {kl/2Ml) D^j,. Summing over aU produced hadrons 
h and integrating over z one can use completeness in the final state to see e.g. that one has (for any flavor) the 
momentum sum rule J dz zDi{z) = 1. 

A note on the names of the fragmentation functions. Because of the analogy with the distribution functions 
in essence the same names are used but using the capital letter, except for the functions corresponding to 
the /.. distribution functions, for which we have used the more familiar name D . The connection with the 
naming scheme used by Jaffe and Ji, namely using the same letter as the distribution functions with a hat is 
straightforward. We have used capital letters in order to avoid double symbols over the function names. 

The twist-two functions have a natural interpretation as decay functions, A^''' 1 being the probability of a 
quark to produce a hadron in a specific spin state (characterized by the spin vector Sh)- Note that the decay 
probability for an unpolarized quark with nonzero transverse momentum, can lead to a transverse polarization 
in the production of spin 1/2 particles. This transverse polarization is orthogonal to the quark momentum. 
Such an effect violates time reversal invariance. For an out-state this is induced by nonvanishing interactions 
in the final state. Indeed, as we will see below, the function D^j, is a purely interaction-dependent function. It 

has no equivalent in the distribution functions. The correlation functions A^^ ''^l and AI**^' are differences 
of probabilities for quarks with different chiralities or transverse spins, respectively, to produce a hadron in a 
specific spin state. For the latter also a new decay function appears because of the non-applicability of time 
reversal invariance. A transversely polarized quark with nonzero transverse momentum can produce unpolarized 
hadrons, in particular it can produce spinless particles, such as pions. The relevant function is the one 
appearing in the single spin Collins effect |2j] . Note that neutral pions cannot be produced, since CP-invariance 
of the strong interactions enforces T-invariance. 

A useful fragmentation function is the one for a quark (with momentum k) producing a quark with momentum 
p {p^ — zk^), for which one has A^'' 1 = i 5{1 — z) S'^{pj, — kj,). This result is useful to find cross sections for 
production of a current jet with (measured) transverse momentum pj, in lepto-production. 

The twist analysis for the projections of the quark-quark-gluon matrix elements. 



Af\z,k'^)= I ^10^e^>^<Tr{O\CiO,C;n+)m^D''iOaWmm 



(66) 

?- =0 



is analogous to the distribution part. The correlation function A^'^' is because of the choice of the path in the 
link operator trivially related to k~ A^^l. Using the equations of motion one obtains the relations 

QTapArj '^-erafSi^A E - m Di - I H + t——H^ 

\ z z Mh 



1 rr I ■ ^ 

, aha' 



r, = i?s-mGis + i Eg, (68) 



(69) 

y ^ / 

Again it is useful to consider the quantity 

A^I^l (z, fe^) = A^Pl (z, - fc" AFI (z, fe^), (70) 

which vanishes for a = — and which for transverse indices reduces after gauge-fixing to the quark-quark-gluon 
correlation function (Fig. 3) needed in the calculation of the Icptoproduction amplitude. As for the distribution 
functions, it is used to define interaction- dependent parts and to express the twist-three functions into twist- two 
fragmentation functions and interaction-dependent parts, explicitly given in Appendix C. 
A consistent definition of antiquark fragmentation functions implies 

A,,{Ph,Sh;k) = -A,jiPh,Sh;-k), (71) 

with the symmetry relations Di{z,k'^) = Di{—z,k'^) and identically for Git, Hit, H^, Et, Gj^, H and Hl, 
while GiL{z,k'^) = -GiL{-z,k'^) and identically for D^j,, H^^, H^, E, D^, El, G'j,, G^, and Ht- 



IV. THE COMPLETE TREE-LEVEL RESULT 

Up to 0(1/(3) one needs to include the contributions of the handbag diagram (Fig. 4), now calculated up 
to this order with in addition irreducible diagrams with one gluon coupling either to the soft part involving 

hadron H or the soft part involving hadron h (sec Fig. 5). The expressions thiis involve the quark-quark-gluon 
correlation functions. The momentum-conserving delta-function at the photon vertex is written (neglecting 
1/Q^ contributions) as 

5'^{p + q-k) = S{p+ + q+) 5{q~ - k') 5^{p^ + Qt - k^), (72) 

fixing p+ = — g+ = XbP^ and = q~ = / zu- The full result, neglecting l/Q^ contributions and still 
omitting flavor indices, is then given by 

-Tr (p)^^A(fc)) - Tr (^^^ A_^^A(fc)7.$^t(^)^ 

-Tr (^7,^7„$(p)7^A^t(fc)) - Tr (^7,^7^AS(fc)7.$(p)) | • (73) 

In this expression the terms with fi^ arise from fermion propagators in the hard part neglecting contributions 
that will appear suppressed by powers of , i.e., 

^1 ^ + m _ {pj + q+Yr 
{pi + q? ^ 2{p+ + q+)q- 

(fci - qf ~ -2(fcj- - q-)q 

Furthermore, the only relevant quark-quark-gluon correlation functions are those in which the argument of the 

gluon field is equal to that of one of the quark fields, discussed in the previous section. These correlation fimctions 
can be related to the quark-quark correlation functions using the QCD equations of motion. Such relations are 



7 

2q- 



+ {-2q+) QV2^ ZhQ^ 



(74) 
(75) 



essential to ensure that at order 1/Q the sum of diagrams in Figs 4 and 5 yields an electromagnetically gauge- 
invariant result. 

The result can be expressed in terms of the twist- two and twist-three distribution and fragmentation functions 
and perpendicular tensors and vectors. These vectors, k'^, p^, 5*^ and /Sl^^, are obtained from the expansion 
into the Cartesian directions of fc^, p^, S!^ and S^rp, respectively. For that one can use the result 



(76) 



obtained from Eqs ^ and The full expression for the symmetric and antisymmetric parts of the hadronic 
tensor are, 
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(78) 



where {/i v} indicates symmetrization of indices and [/i v\ indicates antisymmctrization. Note that one can 
make the replacement f I'^e'^'^a^p e^^'^P'^ ai_pq„. The quark distribution functions in hadron H depend on Xb 
and p^, while the quark fragmentation functions into hadron h depend on Zh and \Phi^ ~ ^PtP, the squared 
perpendicular momentum of hadron h with respect to the quark. 

The cross sections are obtained from the hadronic tensor after contraction with the lepton tensor (Eq. U). 
These contractions, given in Table || use azimuthal angles in the perpendicular plane (see section 2) defined with 
respect to the (lepton) scattering plane and the (spacelike) virtual photon momentum, 



x ■ ax 



'X ■ aj_ 



e^Xf^a^i, = xAa_L = a,^ = |ax| sm0a, 
where we have used a_L A bx = e1'^a±p&±(T- 



(79) 
(80) 



V. RESULTS FOR LEPTOPRODUCTION INTEGRATED OVER TRANSVERSE MOMENTA 



After integration over the transverse momenta {Ph± = —Zh^r) the integrations over k^ and p^ in Eqs 
and can be performed leading to 
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fiXh Gi + XgiDi 
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where the quark distribution functions in hadron H depend on Xb, while the quark fragmentation functions 
into hadron h depend on Zh- The functions indicated with a tilde appear in the quark-quark-gluon correlation 
functions and discussed in section 3. These interaction-dependent correlation functions arc explicitly 
given in Appendix C, just as their fc^-integrated results. The function Df'^^\z) is the (fe^/2M^)-weighted and 
fe^-integrated result of D^j,{z,k'^) (see also Appendix C). 

Separating the cross sections into parts for unpolarized (O) and longitudinally polarized (L) leptons and 
unpolarized (O), longitudinally polarized (L) or transversely polarized (T) hadrons in the initial state one 
obtains the lepto-production cross section 



dxB dydzh 

with (including again the flavor indices) 
daoo 4770;^ s 



da{iH ^ i'hX) daoo _^ daoL _^ daor 



dx B dy dzh dxB dy dzh dx b dy dzh 

^ daLO , daLL , daLT 



dx B dy dzh dxB dy dzh dx b dy dzh 
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The familiar inclusive eross section is easily obtained by using the results for the fragmentation function of a 
quark into a quark and summing over the quark spins, 



da{£H i'X) daoo _^ ddLL _^ daLT 



dxg dy dxs dy dxg dy dxg dy 



with 



daoo Atto^ s 



dxgdy 



(y + 1-2^1 

a, a ^ ^ 



dxgdy 

= ^Xe\S±\ ^e^2y^Jl - y cos{(t>s) g^ixs). 
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(90) 
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These inclusive results lead to the tree-level results for the structure functions Fi{xb,Q'^)-i F2{xb,Q'^), 
g^{xB,Q^) and g^{xB,Q^) which were given in the Introduction. Comparing semi-inclusive and inclusive 
results one finds the familiar result for the number of produced particles 



Nh{xB,Zh) 



dxBdydzh _ X^a.a fli^s) DljZh) 

daoo ~ f^fT \ 

2^a,a^aJl{XB) 



(93) 



Asymmetries in the number of produced particles for polarized leptons and/or polarized target hadrons are 
obtained from the cross sections by putting Sh to zero, 

Nh{xB,Zh;\e;X,S±) = Nh{xB,Zh) [1 + \S±\sin{(l)s) ANh,oT + Xe>^^Nh,LL + Xe\S±\ cos{(j)s) ANh,LT] , (94) 
involving the nonvanishing ones of the following asymmetries. 
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For the double asymmetries, ANh.LL and ANh^lt it is (at least in principle) sufficient to consider either the 
asymmetry for the electron spins or for the target spins since ANh.oL = ANh,LO — 0. Although ANh,OT 7^ 
0, the asymmetries in ANh,OT and ANh,LT involve different transverse directions. Of the asymmetries in 
the number of produced hadrons only AN^^ll is leading. The other asymmetries are proportional to l/Q. 
For a target with transverse polarization in the scattering plane the result is proportional to the transverse 
spin distribution hi and the twist-three fragmentation function H"". For a target with transverse polarization 
orthogonal to the scattering plane the result contains the product of the twist-three distribution with the 
ordinary fragmentation function and the product of the transverse spin distribution hi with the twist-three 
fragmentation function E'^ |l9|. We note that in the naive parton model approach the tilde functions do not 
appear at tree level. 

Measuring polarization in the final state, gives several new possibilities to study quark-quark and quark- 
quark-gluon correlation functions. As discussed in Appendix A the coefficient of Sh in da/dxBdydzh determines 
the polarization. It is convenient to parametrize the polarization vector V^^ in the same way as 5','^' (Eq. |55| ) 
with lightcone helicity A^ and transverse polarization "Phi^ ■ The polarization induced in the final state is 
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In leading order the final state polarization is proportional to the corresponding polarization (longitudinal or 
transverse) of the target hadron and the polarization (longitudinal) of the lepton. The transfer parameters 
involve the leading distribution functions /f , g1 and hi and the leading fragmentation functions and iff. 
For instance, the measurement of transversely polarized A's has been proposed |jl^ as a way to obtain the 
transverse spin distribution hi . 

In the complete result up to order l/Q final state polarization arises from several other polarizations in the 
initial state. Especially noteworthy is the (transverse) polarization in the production of a spin 1/2 particle (e.g. a 



A) for an unpolarized initial state, the term involving /^{xb) Df^^''°'{zii) in Eq. p^02|. This polarization is purely 



transverse, orthogonal to the lepton scattering plane. It involves the (interaction-dependent) time-reversal odd 
fragmentation function D^j?. 



At order 1/Q one finds also transverse polarization in the production of hadrons, coming from longitudinal 
polarization in the initial state or vice versa. This offers new possibilities to obtain twist-three distribution 
functions such as g^{x) and also h1{x). For the latter this semi-inclusive polarization transfer measurement is 
the equivalent of a double-spin asymmetry measurement in Drell-Yan scattering ||l9|,^ . We note the presence 
of several interaction-dependent functions which complicate the naive parton picture for these observables. 



VI. AZIMUTHAL ASYMMETRIES 

The transverse momentum dependence in Eqs|7^ and|7^ lead after integration over and to a dependence 
on Qrp = ~Ph±/z. Furthermore transverse directions appear in the spin vectors. It is useful to project the 
momenta p± and k± onto the direction h, 



{h-a^)h'' + {hAa^)e''/hp, 



(103) 
(104) 



We will not consider polarization in the final state. The result for the symmetric and antisymmetric parts in 
the hadronic tensor are, 
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with (including again flavor indices) convolution integrals of the type 
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For a gaussian transverse momentum dependence the relevant integrals are evaluated in Appendix D. In principle 
one can now identify structure functions by looking at the most general expression for the hadronic tensor pl| . 
We will omit this step. Rather, we immediately give the cross section following from the above tensor. We will 
use the results for gaussian transverse momentum dependence. The complete result for the cross section is again 
written as a sum over terms depending on lepton (O and L) and target hadron (O, L and T) polarizations. The 
contributions that arc nonzero in leading order in l/Q are 
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where G{Qt', R) = {R'^ 1^^) exp{—Q^R^), i.e., a gaussian of which the faU-ofF is determined by a radius R. This 
radius is related to the radii R^ and R^, governing the fall-off of /(.t, fc^) and D{z, k'^) as R^ — R^RfJ {R^+Rf,)- 
These radii again may depend on the longitudinal momentum and on the specific function, i.e., Rh — R^{xb) 
and Rh = R^{zh). Note that in leading order (Jlo = 0. These results have been discussed in Refs p3| and 
psf . The possib ility to measure hi by considering asymmetry in the production of hadrons (e.g. pions), the 
first term in Eq. 110, was first discussed by Collins The complete results show that the measurements of 
specific azimuthal asymmetries in deep-inelastic lepton-hadron scattering with polarized beam and target allow 
the measurement of all six {x- and fcr-dependent) quark distributions, of course keeping in mind that also a 
separation of the different flavors is needed. For obtaining a first indication for the behavior of the results, the 
dominance of u-quarks in the proton, especially at moderate a;- values will be helpful. In all the expressions only 
two fragmentation functions play a role. The first is the 'ordinary' fragmentation function Z?", the second is 
the time- reversal odd function ifj*-" . 

In principle the general expression for the hadronic tensor in the beginning of this section allows the calculation 
of all azimuthal asymmetries in the cross section that contribute up to order 1/Q. Below, we give first the results 
for unpolarized spin 1/2 targets, which for unpolarized leptons is the extension of Eq. 108, 



daoo 



dxs dy dzh d'^P 



-2(2-y)yT^cos(</),)^ 




■xlf^-{x,)Dl{zn) 



-^-b/i (-^b) 



D^''{zh)\ Q{Qt;R) 



Zh 



(113) 



da 



LO 



dxB dy dzh d^Phi_ 



AtTO^S ^ 2o /T ■ A Qt MR'^ 2. MT±a/ .G{Qt\R) ,,,,, 

Ae }_^ei2yy/l - y sm (l)h — ^^j-^ xi e^ixs) H^'^izh) — (114) 



Q MhRi 



The {cos{4>h)) asymmetry in unpolarized leptoproduction, unfortunately is rather complicated, involving one 
twist-three distribution function (/^°) and one twist-three fragmentation function {D^°-) |jl^. It is important 
to point out, however, that the (cos(0/i)) asymmetry is not only a kinematical effect. It reduces to a kinematical 
factor only depending on y and when the interaction-dependent pieces in the twist-three functions are set 
to zero, /^'^ = and D^" = 0, implying f-^" = f^/xs and D-^" = Zh D^f- At order 1/Q there is no {cos(2(j}h)) 
asymmetry in the deep-inelastic leptoproduction cross section. For polarized leptons and unpolarized targets a 
(sin((/)/j)) asymmetry is found |^^, involving the interaction-dependent part of the distribution function e° and 
the time- reversal odd fragmentation function H^'^. Noteworthy is that it is the same fragmentation function 
that appears in several of the leading azimuthal asymmetries for polarized targets. 

As indicated, the above asymmetries have been discussed in various papers. The remaining asymmetries at 
order 1/Q, i.e. those for polarized targets, can also be extracted from the general result at the beginning of this 
section. They are 



dtJoL 



dxg dy dzh d'^Phi_ 



— — \Y,elU\-y) sm(20,) 



MMh Rl Rl 



Xshi^{x^)Ht%Zh) 



-2(2 -y)^!^ sin((/),0 
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Q [MMhR^Rt 
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xlhl{x^)Ht%Zh)- 



MhR^ 
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H ^ J \ 



a. a 



+(1 - y) sin(30k - 4,,) 2M2MXi Rl ''•''■^T{''.Wt'M 



QlR^ ^ .H^\ 1 g(QT;i?) , . 



+2(2 - 2/)yr^ sin(20, " '/'^) ( ^ 2Ah Rj,Rl l^'^i^^-) ~ a:,/i^/(.T,) j H^'^iz^) 



a, a 



MhR"" ^^a, .E''{zh)\\G{QT\R) 



^ x^h^{i(x^)—-^ ) ) ^^^^i:^, (117) 



dxsdydzhd^Ph± fl' " r ^ 2/ ' ^ ' MR]j 



+ 



[l-Q|i?2] x«5°^(xb) 



Q 2MMhR]jR:i ^ ) Zh 



Qt QtR' 2 ±ar ^^a, ^, QtR^ Mh ,E-izh) 



-2yy/l^ cos(20;, " "^^^ ^ I ^-5t ^?(^^.) + ^ xMri^ 



Zh 



2MRIR1^^3^- Zh )] zl • 

Because the results involve mostly higher harmonics in the 0/i-dependence and involve most of the allowed 
combinations of twist-two and twist-three distribution and fragmentation functions discussed in section 3, we 
do not expect that these results soon will be used to extract the new information they contain on the structure 
of the nucleon, specifically the quark-quark-gluon correlations contained in the interaction dependent pieces. 
In first instance, the result have been given because they can be used to estimate azimuthal asymmetries. 
Such estimates may be necessary in the absence of full azimuthal coverage. We note that the estimate of the 
asymmetries in the naive parton model is obtained by setting all interaction-dependent pieces (functions with 
a tilde) to zero. In that case the twist-three functions can be expressed in the twist-two ones. The explicit 
relations are given in Appendix C. Note, however, that although the functions appear in the 1/Q parts of the 
cross section, there is no reason to expect that the interaction-dependent functions themselves are smaller than 
the twist-two functions. They are just different matrix elements involving gluon fields. 

Finally, we want to discuss the case in which only the azimuthal direction of the current jet is detected, 
without an analysis of this jet. The measurement of the (small) transverse momentum of the current jet allows 
(in principle) a direct study of the transverse momentum of the distribution functions. The starting point are 
Eqs 13 and [7q in which one uses the quark — > quark fragmentation function. The nonvanishing contributions 
in the cross section (with p_L being the jet transverse momentum involving azimuthal angle 0^ ) are 



dxs dy d^pA 



d(JLL 47ra^ s 



dxs dy d^pA 
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-2(2-y)v/T^cos(^,)^4/^n^B,Pi)|, 
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(119) 



(120) 



a, a I, 



'2yy/l - y cos(0s) ^ xig^{xs,p\) 



-2yy/T^ cos(2</., - 0,) — |^xig^na;.„pi)|. (121) 

Again decomposing the twist-three functions in these jet asymmetries into their twist-two parts and the 
interaction-dependent piece, one obtains the naive parton model results, if the interaction-dependent pieces 
are set to zero. For the unpolarized case, the cross section then becomes proportional to /f , reducing the 
(cos(^j)) asymmetry to a purely kinematical effect depending only on y and Q as discussed in Ref. Q. Simi- 
larly, in the case of longitudinal polarization setting the interaction-dependent pieces to zero, produces a result 
proportional to the quark helicity distribution gi plus a quark mass term. Except for this mass term the 
(cos((/)j)) asymmetry is purely kinematical in that case. We would like to turn things around, however, and 
state the importance of finding deviations of the naive parton results for unpolarized and polarized jet pro- 
duction in order to get an idea of the quark-quark-gluon correlation functions f^"^ and gj^"^ in a nucleon. For 
transversely polarized hadrons one can study g!^(xB,p'^), which integrated over p± can also be measured in 
inclusive leptoproduction. 



VII. SUMMARY 



In this paper the complete tree-level result up to order 1/Q for polarized deep-inelastic leptoproduction has 
been presented. The formalism is the diagrammatic approach in which soft hadronic parts are represented by 
expectation values of nonlocal - here bilocal - combinations of quark and gluon fields, referred to as correlation 
functions. Both quark-quark and quark-quark-gluon correlation functions need to be included in order l/Q. 
The quark-quark-gluon correlation functions can be related to the quark-quark correlation functions using the 
QCD equations of motion. Such relations are essential to ensure an electromagnetically gauge-invariant result. 
Essential in our treatment is also the explicit treatment of quark transverse momenta in order to study azimuthal 
asymmetries in semi-inclusive lepto-production processes. 

The full result for the hadronic tensor contains a large number of terms. They are some of the allowed terms 
in the most general expansion of the hadronic tensor, which can be expressed in terms of 'standard' tensor 
structures built from the momenta and spin vectors of the hadrons involved multiplied with structure functions. 
Having such an expansion one could read off the result for the structure functions. In view of the large number 
of structure functions that are allowed in principle, we found it more practical to split up the cross section in 
the parts involving the lepton polarizations (unpolarized or longitudinal polarization) and hadron polarizations 
(unpolarized, longitudinal and transverse polarization). 

In the leptoproduction cross section in which one integrates over all transverse momenta the results have been 
presented as the numbers of produced hadrons Nfi{xB,Zfi) and the asymmetries in these numbers arising from 
the polarization of initial state lepton or hadron. Finally the polarization of a produced (spin 1/2) hadron has 
been discussed. The results include among others the following features 

• An asymmetry in produced particles (e.g. pions) is obtained at order l/Q for unpolarized leptons and 
transversely polarized hadrons (Eq. ^6|) proportional to the transverse spin distribution and a twist- 
three fragmentation function. 



A transverse polarization orthogonal to the scattering plane is induced in spin 1/2 hadrons in the final 
state (e.g. A-baryons) starting with an unpolarized lepton and unpolarized target hadron. This polar- 
ization is proportional to the unpolarized quark distributions /f and the fragmentation function D'^^^ " 



(Eq. 102). For longitudinally polarized lepton and target hadron a similar polarization emerges, now pro- 



portional to the quark helicity distribution g° and the same fragmentation function D^^^ °. Although the 
fragmentation function is a new one, allowed because of the non-applicability of time-reversal invariance 
in the production, the comparison of induced polarization for polarized and unpolarized initial state may 
give another handle on determination of quark helicity distributions, especially because the production of 
A's produces different flavor weighting compared with inclusive processes. 

For completeness also all leading as ymmetr ies in number of particles and in induced polarization have 



been included in Eqs |95|-|99| and Eqs 100-102. These are the ones that appear in the naive parton model. 



The dependence on quark transverse momenta becomes explicit in semi-inclusive lepto-production in the 
transverse momentum of the produced hadrons. Here the following features emerge, 

• All six twist- two X- and p^^-dependent quark distribution functions for a spin 1/2 hadron can be accessed 
in leading order asymmetries if one considers lepton and hadron polarizations. One of the asymmetries 
involves the transverse spin distribution On the production side, only two different fragmentation 
functions are involved, the familiar unpolarized fragmentation function Df and the (interaction-dependent) 
fragmentation function _ffj^°. 

• At order 1/Q a (sin(0^)) asymmetry is found for polarized leptons and unpolarized hadrons. This probes 
the interaction-dependent distribution functions e° in combination with the same fragmentation function 
Hi'^'' mentioned in the previous item. 

Several more azimuthal asymmetries are found for polarized targets. The results have been given explicitly 
because they indicate the type of asymmetries to be expected in deep-inelastic leptoproduction. At this point 
we also want to point out once more that our results represent tree-level calculations. The inclusion of gluon 
ladder-graphs give corrections proportional to as and as InQ^. For the cross sections integrated over trans- 
verse momenta, they lead to scale-dependence of the distribution functions. They will also enter as additional 
contributions to several of the observables, including the ones that are zero at tree-level. A well-known example 
of the latter is the longitudinal structure function Fl ~ F2 — 2xbFi, which is proportional to a^. Work along 
these lines with explicit treatment of transverse momenta of quarks is under investigation. 

We acknowledge contributions of J. Levelt (University of Niirnberg-Erlangen) in the early stage of this work. 
We thank D. Boer and R. Jacob for a careful reading of the manuscript. This work is part of the research 
program of the foundation for Fundamental Research of Matter (FOM) and the National Organization for 
Scientific Research (NWO). 



APPENDIX A: SPIN VECTORS 



In this appendix we would like to remind the reader of the role of the spin vectors S and Sh for initial state 
and final state hadrons, respectively. A more extended discussion of aspects related to the role of the spin vector 
of the produced hadron (applied in e^e~ annihilation) can be found in Ref. ]2^ . Starting with the quantities 

^t^M^P^Ph) - {h;k\T\H;t){H;j\T^h;l) 

— generalized tensor in 2 x 2 H-spin-space (ij) and in 2 x 2 h-spin-space (kl), 

/r, i-TT -iV^i X / x; x H rest frame 1/, „\ 

Prj{P,S) = {H-i\Y,\a)p{a){a\H-j) = -[l,, + a,, ■ Sj 

a 

— initial state density matrix , defining 5, 
RhiP.J) = {h;k\f){f\h-l) ^'^e^U'-^™^ w{f)(l,i+alrS,) 

— decay matrix into final states f in 2 x 2 h-spin-space (kl), defining Sh, 

where the latter is normalized to X]/-^fe;(/) ~ ^ki- Note that the summation over / is just symbolic for all 
kinematic variables appearing in the decay. An example of w{f) is the center-of-mass distribution w(6'^Ar) 



in the decay of the A baryon. Starting with some (general) initial state, one has the following result for the 
semi-inclusive hadronic tensor in which the decay products of hadron h are detected, 

- /) = PiiiP, S) P, Ph) RikiPh, f) = P^ S; Ph) RikiPh, /). (Al) 

This shows for instance why the spin vector, originating from p{P, S) can only appear linearly in the parametriza- 
tion of the tensor. Summing (or integrating) over the decay products of h one is not able to measure the 
polarization of h, 

Y.^^'i.i^h^ f)= W^'^{q; P, S; Ph) = W^-{q; P, S; Ph,SH = 0). (A2) 
/ 

The general result in the h-spin-space can (in the rest-frame of hadron h) be written as the unpolarized result 
times the production density matrix parametrized in terms of the polarization P^, 

Wfc7(g; P, S; Pu) = W^'iq; P, S; Ph, 0) i (l;^; + aui ■ Vh{q; P, S; Ph)) ■ (A3) 
Thus one can write in the h-iest frame, 

yV^-'ii - /) = W^^iq; P, S; Ph, 0) w{f) (l + r{q; P, 5; Ph) ■ Sh{Ph; /)) 

= W''''{q;P,S;Ph,Sh)w{f). (A4) 

This shows that also the spin vector 5*^ appears only linearly, while it multiplies the polarization vector as 
appearing in the production matrix. Parametrizing the polarization vector Vh similar as Sh in a longitudinal 
polarization Ah and a transverse polarization VhT, which because the polarization vector parametrizes a density 
matrix satisfy A^ + < 1, one can in a general frame write for V ■ Sh = AhXh + 'Ph± ■ Sh±- 



APPENDIX B: THE LINK OPERATOR 



In this appendix we discuss the path structure of the link operator used in the definition of the various 
correlation functions. In hard processes one encounters the correlation functions integrated over at least one 
lightcone component of the parton momentum. The actual component which one needs is actually fixed by the 
large vector q. This vector can be used to fix a lightlike vector for each given (small, i.e. <C Q^) vector p, 

For the correlation function $(P, S;p) one obtains as a consequence of the fact that p^ ^ P^ ^ p ■ P <^ 

np^np= 0, Ot) = n_. (B2) 

The lightlike vector is (approximately, i.e., up to corrections of 0{1/Q'^)) equal for all momenta in a 'connected' 
soft part. These momenta themselves define the 'other' lightlike vector, p/Q oc P/Q oc n+. Similarly one finds 
for A{Ph,Sh;k) that 

nfc«np, = [o, ^, Ot) = n+, (B3) 

while k/Q (X Ph/Q oc ri-. The lightlike vector determines the integration used, for instance in J dp~ $(P, S\p) 
the component p~ = {Q/A)p- n+. It also fixes the path structure in the link operator 

C{0,x;n-)=Pexp(^-ig dsf" A^{s)'j , (B4) 

The link operator £(0,x;n_) is actually constructed as the average of two links with paths Pi and P2. For i> 

both lie in the a;"^ = plane and run along the direction except for endpoints at x~ = ±00. The path Pi runs 
from [0,0, Ot] to [x^,0,Xt] via [-|-oo,0,Ot] and [+oo,0,Xt], while the path Pj runs from [0,0, Ot] to [x^ ,0,Xt] 



via [—00, 0, Ot] and [—00, 0, x^]- In the gauge y4.+ = ?i_ • A = 0, the link operator in the projections $^1 becomes 
unity, except for the parts at x~ = ±00. However, it is well-known that = does not completely fix the 
gauge. The residual gauge freedom is for instance fixed by imposing an antisymmetric boundary condition on 
At (HJ. Fixing in this way the residual gauge freedom on At, ^t[— 00, cc^] = — ylT[oo, a;^], allows a 
unique inversion of the relation G"*"' = 9"^^^, 

J^T\^-,x^,x^\ = \^dy-t{x- -y-)G^\y-,x^,x^\ = ^G+\x- ,x-^ ,x^\. (B5) 

The choice £(0,a;;n_) = (£i(0, a::;Pi) + £2(0, a:; P2))/2 then reduces to unity. 

For the link operator in the correlation functions A all choices are based on the n_(--direction. It becomes 
simply unity in the gauge fixed by A^ — and an antisymmetric boundary condition for At- 



APPENDIX C: EXPLICIT QUARK-QUARK-GLUON CORRELATION FUNCTIONS 

For the correlation function. 



Pt) - - '^^^^ Pt). (CI) 



I rpj 

one has $^ = while the functions with transverse indices reduce after gauge fixing (A+ = 0) to a quark- 
quark-gluon matrix element (Eq. |5^). It is convenient to parametrize the following combinations (in which all 
indices a and /3 are transverse only) as 

3™/3$f"'^^^'=A/x/i„ (C3) 
$^[7+1 _ ,,^/3 $ j^^+^^l ^ - ie« x~g^ - zefSTp Mx cj't, (C4) 

in terms of (interaction-dependent) functions indicated with a tilde. Using the results for as following from 
the equations of motion (see section 3), one can split the twist-three distribution functions into the twist-two 
part and an interaction-dependent part. The first part is the one which is obtained in the naive parton model. 
The complete list is 

e{x,pl)^^ll^^ + e[x,pll (C5) 
M X 

f^(x,pl) = ^-^^^+~f^{x,pl), (C6) 

^ rn hMx,pl) ^ ^^^^ 
M X 

giix,pl) = ^l^i^ -h ^ + giix^pll (C8) 

X M X 

g-ix,pl) = ^-^I^ + ^ +^^(.,p^), (C9) 

h^ix,pl) = + h^ix,pl), (CIO) 

, / 2 ^ m giL{x,pl) pI hj;^{x,pl) .7,9. ,^-1 -, s 

hL{x,p^)^— — + hL{x,p^), (Cll) 

. / 9. m giTix,pl) hiTix,pl) pI h^j,{x,pl) ~ ^ 

hT{x,p,) = ^ + hT{x,p^). (C12) 

The Pr-integrated results involve the functions 

e{x) ^ J d'pr e{x, p2 ) ^ ™ ZlM + g(a;), (C13) 



9t{x) 



9rri^) , m hi{x) 



M X 



-9t{x), 



M X 



(C14) 
(C15) 



where the functions e{x), gxix) and h^^x) are the pT-integrated results of the interaction-dependent parts. The 
functions with upper index (1) indicate p^/2M^-weighted functions, 



/2 

ht^'\x)^ Jd'p,^hi,{x,pl). 



(C16) 
(C17) 



Using the exphcit expansion into amphtudes one can derive the foUowing relations for the px integrated distri- 
bution functions pC[|, which clearly exhibit the role of transverse momenta and interaction terms, 



grix) 



^ , qi(v) m 
V M 

9i{x) + ^5(t(2^) 



hi{x) 



\ hi{y) 



■9t{x) 



hL{x) = 2x I dy 



hijy) rn 
f M 

h,{x)-2-^-ht^'\x). 



3M-2X tdy'-^ 



y" 



+ hLix)-2x i'dy^ 



(C18) 



(C19) 



The decomposition of gx contains the Wandzura- Wilczek part p9[ , a quark mass term and interaction-dependent 
parts. Noteworthy is the appearance of the (p^/2M^)-weighted twist-two functions, e.g. the function g[l^(x) 
[In Ref. referred to as B^{x)]. 
For the fragmentation part one has for the projections 



A^[^l (z, k^) ^ (z, fc,) ~ fc" API (z, fc,) 



the result A = 0, while for transverse indices 



4z (27r)3 



Tr (Om) A"(x)ala,^(0)r|0) 



(C20) 



(C21) 



It is convenient to parametrize the following combinations (in which all indices a and P are transverse only) as 



.Plia" 75] 



Z Z 



,a/3 A [7 75] _ .a 



M/3 



1 ^ 



(C22) 
(C23) 

(C24) 



in terms of (interaction-dependent) functions indicated with a tilde and depending on z and fc^ = —zk^.- One 
again can use them to split the fragmentation functions into a 'parton' piece and interaction dependent functions. 



D^j'{z,k'^) — D^rp{z,k'^), 



(C25) 
(C26) 



E{z,k'' 



D^{z, k'^) = zDi{z, k'^) + D^{z, k'^), 
E,{z,k'^)=Es{z,k'^), 

G'Az, k'^) = ^ zHMz, k'^) + G't{z, k'^), 
G'^{z, k'^) = zGis{z, k'^) + — — zH^g{z, k'^) + G^(z, fc^), 
H^{z,k'^) = zHiT{z,k'^) + H^{z,k'^), 



H{z,k'^) 
Hs{z,k'^) 



^^zHtiz,k'^) + H{z,k% 



m 

ii4 



zGis{z, fc'y) - 



''^ ■ ^''^ zH,t{z, k'^) - ^ zHtiz, k'^) + Hs{z, k'^) 



Ml 



(C27) 

(C28) 
(C29) 
(C30) 

(C31) 
(C32) 
(C33) 

(C34) 



Note that the twist-two time-reversal odd functions D^j, and are interaction-dependent. Below we give 
some of the fe^-integrated functions that turn out to be relevant in the integrated cross sections. They are 



E{z) = / (fk'^ E{z, k'^) = z^ \ (fkj. E{z, z^kl) = zDi{z) + E{z), 



Gt{z) = / (fk'^ 



G't{z, fey ) + 2^ Grp{z, k'^) 



Mh 

= zG[%) + ^zH,{z) + Gt{z), 



Hl{z) = I d'k'^HL{z,k'^) = ^2zH^l'\z) + --zGiiz) + HLiz), 
J Mh 

while the fe^/2Af,j-weighted functions (indicated with index (1)) are, 



Dt^iz) 



G^l^iz) = I d'k 



j d^k'^ 2M^ ^it(^' '''t)' 



2Ml 



GMz.k'^), 



(C35) 
(C36) 
(C37) 

(C38) 
(C39) 
(C40) 



The above decomposition leads in combination with the amplitude expansion (based on Lorentz invariance) to 
the following relations, 



Gt[z) = z I dy —2 h — 



Mh 

Gi(z)-Hz^(zG1^3^(z)), 



zHi{z)-z l\y^^^^^ 



y 



+ Gt{z)-z l\y^^^'"^ 



1 ' 



y Mh 
d 



zGi(z)-2 / dyGiiy) 
H,{z)^2z^^(zHtl^\z 



HUz)-2l\y^ 

y 



(C41) 



(C42) 



Analogous to the sum rules dx (72 (a;) = | |30| , where (72 = 5t — ffi , and dx /12 (x) = ]3l|j20| ] , where 
h2 = 2{h]^ — hi), one obtains sum rules 



where G2 — Gt — Gi and 



(C43) 



(C44) 



where H2 = 2{Hl — Hi). As for the distribution functions, one must be aware of convergence problems in the 

integrals and of the fact that in the cross sections other contributions beyond tree-level may be present. We 
consider these sum rules at this point as merely academic. 



APPENDIX D: CONVOLUTIONS FOR GAUSSIAN DISTRIBUTIONS 



In order to study the behavior of the convolutions of distribution and fragmentation functions it is useful to 
consider gaussian distributions, 



f{x,pl) = fix,0) cM-RIpI) 
D{z,k'^) = D{z,Q) expi-Rlkl) 

TT Z 



D{z) eM-Rlkl) = ^ g{\kT\;Rh) = D{z) Q ( z\k-r\ 



Rh 
Z 



(Dl) 
(D2) 



In that case the convolution becomes 



I[fD]{x,z,QT)= J cPp^f{x,pl)D{z,\Ph±-zp^\'') 

fix,OT)D{z,OT) 



R'h + RI 

= f{x)D{z) 



R 

Q{Qt;R) 



(D3) 



where B? = Rfj R\/ {Rfj + Rf)- The other convolutions that appear in the cross sections are of the form 
h - Pt 
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h ■ kr. 
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fD 
fD 



MMh 
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M2 

2{h-k.rf-kl 
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fD 



Qx R'^ 
M2 Rjj 



nfD], 



I[fD], 



2{h ■ Pt) {h ■ kT) — Pt ■ kT 



MMh 



fD 



MMhRjjRl 



I[fD], 



2{h ■ Pt) jh ■ fc-r) -pT -kT pI 

MMh M2 

A{h ■ pT)^{h ■ kT) - 2{h ■ Pt) {pt ■ kT) -p\{h- kT) 



O d2 d2 



M^M, 



h 



M^MhR%Rl V R"^ 



fD 



R^ 



MmhR% Rl 



Q'tRIJ I[fD], 
I[fD]. 



(D4) 
(D5) 
(D6) 
(D7) 
(D8) 
(D9) 
(DIO) 
(Dll) 
(D12) 



Note that we have suppressed possible dependence of the radii on the specific function or the kinematic variables 
irrelevant in the integration over transverse momenta. We can trivially generalize the results by taking into 
account such dependence, Rh — * Rh{x) and Ru — > Rh{z). 
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TABLE I. Contractions of the lepton tensor L^v with tensor structures appearing in the hadron tensor. 



= al''elHxp-(el''axp6x.)flr 



i^(l-j/)|ax||6x| cos(0a + 06) 



(1 - y) \a±\ |bx| sm{(j)a + (l>b) 



a±p 



-^(2-2/)yr^|ax| cosc^a 
^(2-2/)Vr^|ax| sin0„ 
Ae^J/(l-i) 

^2/(1-1) |ox||6x|sin(0,-0„) 
-'^e ^ yVl - sin 0a 

-'^e ^ - COSflia 



FIG. 1. Kinematics for iH scattering. 




FIG. 3. Quark-quark-gluon correlation functions $a (left) and (right) contributing in hard scattering processes 
at subleading order. 




FIG. 5. Quark diagrams contributing to IH scattering at order 1/Q. 



